% X0=[20;80;390;0.52;0.41;0.1];%给定初值，不能乱给，会迭代发散，求解失败
% esp=0.0001;%精度
% E=1;%差值
% while abs(E)>esp
%     f=double(subs(fun(X0),{'X','Y','Z','A','B','G'},{X0(1,1),X0(2,1),X0(3,1),X0(4,1),X0(5,1),X0(6,1)}));%subs进行变量替换，将符号变量替换相应数值
%     jacobi=double(subs(dfun(X0),{'X','Y','Z','A','B','G'},{X0(1,1),X0(2,1),X0(3,1),X0(4,1),X0(5,1),X0(6,1)}));
%     X=X0-jacobi^-1*f;
%     E=X-X0;
%     X0=X
% end

fun()
dfun()